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Abstract 

We study the phenomenon of real space condensation in the steady 
state of a class of mass transport models where the steady state fac- 
torises. The grand canonical ensemble may be used to derive the 
criterion for the occurrence of a condensation transition but does not 
shed light on the nature of the condensate. Here, within the canonical 
ensemble, we analyse the condensation transition and the structure 
of the condensate, determining the precise shape and the size of the 
condensate in the condensed phase. We find two distinct condensate 
regimes: one where the condensate is gaussian distributed and the par- 
ticle number fluctuations scale normally as L 1 / 2 where L is the system 
size, and a second regime where the particle number fluctuations be- 
come anomalously large and the condensate peak is non-gaussian. Our 
results are asymptotically exact and can also be interpreted within the 
framework of sums of random variables. We further analyse two ad- 
ditional cases: one where the condensation transition is somewhat 
different from the usual second order phase transition and one where 
there is no true condensation transition but instead a pseudoconden- 
sate appears at superextensive densities. 

PACS numbers: 05.40.-a, 02.50.Ey, 64.60.-i 
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1 Introduction 



Mass transport models encompass a large class of systems wherein 'mass', 
or some conserved quantity, is transferred stochastically from site to site of 
a lattice. Well known examples of such models are the Zero-Range Process 
(ZRP) P, and the Asymmetric Random Average Process (ARAP) |2j. These 
are simple fundamental models which have been used to describe such diverse 
physical situations as traffic flow |3j, clustering of buses [I], phase separation 
dynamics [H], force propagation through granular media [Oj, shaken granular 
gases [7J |Hj and sandpile dynamics [S] . For a recent review of the ZRP and 
related models and applications see |TU] . 

Condensation is manifested in such models when, in the steady state, 
a finite fraction of the total mass condenses onto a single lattice site. At 
low global mass densities the system is in the fluid phase where the mass is 
distributed evenly over all sites. Here, the single site mass distribution p(m) 
typically decays exponentially with m implying a finite amount of mass at 
each site with the mean being the global mass density p. On increasing the 
global density above a critical value, p c , an extra piece of p(m) emerges which 
represents the condensate i.e. a single site which contains the global excess 
mass (p — p c )L where L is the system size. Thus in the condensed phase the 
condensate coexists with the background fluid. Note that the condensation 
occurs in real space and, in principle, in any number of dimensions. In this 
paper we will analyse in detail the structure of the condensate, elucidating 
when it occurs, its shape and its fluctuations. A short communication of 
some of our results has been given in |IT|. 

In general mass transfer models form examples of systems exhibiting non- 
trivial nonequilibrium steady states. Since the models are defined by stochas- 
tic dynamical rules, their steady state distributions are not a priori known. 
However it turns out that a large class of models, which we shall discuss be- 
low, enjoy the convenient property of having a factorised steady state. This 
means that the steady state probability P({mi}) of finding the system with 
mass mi & t she 1, mass m 2 at site 2 etc is given by a product of (scalar) 
factors f{rrii) — one factor for each site of the system — i.e. 

P( W) = Z(M, L)- 1 H fimi) 5 ( mi - M ) . (1) 

i=i \i=i / 

where Z(M, L) is a normalisation which ensures that the integral of the 
probability distribution over all configurations containing total mass M is 
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unity, hence 



Z(M,L) = H / dm,/(m,) * £ 



1=1 L-^o J \ J=1 



mi- M 



(2) 



Here, the 5-function has been introduced to guarantee that we only include 
those configurations containing mass M in the integral. The single-site 
weights, /(to) are determined by the details of the mass transfer rules. 

Thus when the steady factorises, the analysis of condensation reduces to 
the evaluation of (J2J), a problem first addressed in [T2] . Equation (J2J) defines a 
canonical partition function in that the delta function imposes the constraint 
of fixed total mass M. As we shall discuss in detail below analysis of (J2J) has 
only been carried out in the fluid phase and effectively in the grand canonical 
ensemble where the total mass is allowed to fluctuate. Our aim here is to 
present a full analysis of (I1I2|1 in the canonical ensemble. In particular this 
shall elucidate the mechanism of condensation within factorised steady state. 

Although our focus in this paper is on the analysis of the steady state 
(jl!2|) . it is relevant to briefly review some particular models and dynamics 
which give rise to such steady states, and previous studies of condensation 
phenomenon. Firstly we mention the backgammon model ^Hj where unit 
masses hop under dynamics respecting detailed balance with respect to an 
energy function which is simply minus the number of unoccupied sites in the 
system. When the temperature T — > a condensate, where all masses are 
on the same site, dominates the steady state. Since the model and dynamics 
are constructed so that an energy functions exists, the steady state automat- 
ically factorises into the form (II 1211 (but with discrete mass variables). The 
motivation for the Backgammon Model was actually to study a simple model 
for glassy dynamics since in the late time regime at low T entropic barriers 
and slow dynamics arise. 

In ^2] the factorised steady state of the Backgammon model was gen- 
eralised to the form (jip2|) with single-site weights /(to) that could have an 
asymptotic power law dependence 



It was shown that a condensation transition would occur if 7 > 2. Mean- field 
dynamics for this general factorised steady state were constructed from the 
detailed balance condition and the mean-field dynamics of condensation was 
studied in jTJj. Here 'mean field' is used in the sense that a particle can hop 
from a given site to any other. 

On the other hand it has long been known that the ZRP, a stochastic 
model defined by local stochastic dynamics and without respect to any energy 
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function, has a factorised steady state Oil]- I n a heterogeneous system, 
where the rules for mass transfer depend on the site, the condensation may 
occur at the site with lowest outgoing mass transfer rate ^Hl E±] • In this 
case the mechanism for condensation (in real space) is exactly analogous to 
Bose-Einstein condensation in momentum space in an ideal Bose gas and the 
slowest site plays the role of the ground state in the quantum system. 

For homogeneous dynamical rules for the transfer of mass, which we will 
focus on in this work, the condensation mechanism in the ZRP was first 
studied, to our knowledge, in |lj where the ZRP was used as an approximate 
description of a model of bus routes. It was shown that if the hop rate u(m) 
for a mass to move from a site with m masses to its neighbouring site decays 
as u{m) ~ + j/m) then the single-site weight decays as © so that a 
condensation transition occurs if 7 > 2. Condensation in the ZRP has been 
further studied in |THl EES HH EI] • 

Another model with a factorised steasy state is the ARAP [2j. Although 
usually this does model does not exhibit condensation, condensation may be 
induced by imposing a maximum threshold on the amount of mass that may 
be transferred between sites |22j. 

It turns out that the ZRP and ARAP may be unified by considering a 
very general class of one-dimensional mass transport models [23J: a mass m; 
resides at each site i; at each time step, a portion, fhi < rrii, chosen from 
a distribution (p(jh\m), is chipped off to site i + 1. Choosing the chipping 
kernel ip(m\m) appropriately recovers the ZRP, the ARAP and the chipping 
model of [21]. Moreover the class of models encompasses both discrete and 
continuous time dynamics and discrete and continuous mass. A necessary 
and sufficient condition for a factorised steady state is that the chipping 
kernel is of the form 23J 



where u(z) and v(z) are arbitrary non-negative functions. Then the single 
site weight is given by 



This result may be generalised to arbitrary lattices in all dimensions [2*H| EE] • 
Given a chipping kernel (p{fh\m), sometimes it is hard to verify explicitly 
that it is of the form (jSJ) and thereby to identify the functions u(m) and 
v(m) in order to construct the weight f(m). This problem was circumvented 
by devising a test [27 \ to check if a given explicit (p(rh\m) satisfies the con- 
dition (JH]) or not. Further, if it "passes this test," the weight f{m) can be 
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found explicitly by a simple quadrature [27|. Finally, for any desired function 
f(m), one can construct dynamical rules (i.e., </?'s) that will yield f(m) in a 
factorised steady state. 

Condensation has also been noted and studied in mass transport models 
without factorised steady states. In the chipping model defined in [21] all 
the mass at a site can move to a neighbouring site. Condensation in the 
model was studied in a mean-field approximation for the steady state which 
amounts to approximating the true (unknown) steady state by the factorised 
form (jl|2j) . It turned out that in the true steady state condensation is actu- 
ally only observed in the case of symmetric hopping flQ EBl EUJ • However a 
generalisation which includes both the chipping model and ZRP as special 
cases does appear to exhibit condensation for asymmetric hopping [30] . Fi- 
nally we mention more complicated condensation scenarios such as the two 
species ZRP [HI] where, for example, one species may induce condensation 
in the other. 

So far we have reviewed models exhibiting condensation which comprise 
hopping particles. Condensation phenomena have also been observed in the 
context of the rewiring dynamics of networks with a fixed number of nodes 
and links [H21EH]- There the condensate corresponds to a node which is 'hub' 
which captures a finite fraction of the links. Also, models of macroeconomies 
may exhibit 'wealth condensation' where the wealth is the conserved quan- 
tity and condensation corresponds to a single individual capturing a finite 
fraction of the wealthfHI]. A further important application of the conden- 
sation mechanism is to determine a criterion for phase separation in one- 
dimensional systems Phase separation, or indeed the lack of it, has 
been reported in a number of one-dimensional and "quasi- 1-D" systems typi- 
cally involving three species of charged particles with an exclusion interaction 
|Hni EH EE1 EH EH El H2j- By mapping the domains between the neutral 
particles onto the sites of a ZRP, the ZRP can be used as an effective descrip- 
tion of domain dynamics and the criterion for condensation in this effective 
description determines when phase separation will or will not occur. 

The paper is structured as follows. In section El we review the grand 
canonical approach which correctly predicts the critical density, but does 
not shed light on the nature of the condensate. Before embarking on the 
canonical analysis we interpret in section El the canonical partition function 
in the framework of sums of random variables and deduce rather quickly some 
key properties of the condensate. In sectionElwe summarise the results of the 
full canonical analysis which are detailed in sections 03 and El in section El we 
present an exactly solvable case where a closed expression can be found for the 
canonical partition function and in section El the asymptotic analysis of the 
canonical partition function in the general case is presented. In sectional we 
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consider other forms of f(m) such as Q with 7 < 2 for which condensation 
does not occur, and stretched exponential f(m) for which condensation does 
occur but with some differences. We conclude with a discussion in section |H1 



2 Criterion for Condensation: Grand Canon- 
ical Approach 

Assuming that the steady state of mass transport model is of the factorised 
form fjl|2D . one next turns to the issue of condensation within this factorised 
steady state. In particular, we ask: (i) when does a condensation transition 
occur, i.e., the criterion for condensation (ii) if condensation occurs, what 
is the precise nature of the condensate? The factorization property allows 
(i) to be addressed rather easily within a grand canonical ensemble (GCE) 
framework - an approach similar to the one used in traditional Bose-Einstein 
condensation in an ideal gas of bosons, though there are some important 
differences in the present case. In this approach one takes the thermodynaimc 
limit (L,M — > 00 with fixed overall density p = M/L). Then becomes 
a product measure state where the sites are decoupled from each other and 
the single site mass distribution is simply given by p(m) = f{m)e~^ m where 
p is the negative of the chemical potential and is chosen to fix the density, 
giving rise to the relation 



The criterion for condensation can be derived easily by analysing the function 
p(p) defined in Eq. ©. 

The behavior of p(p) depends on f(m). We consider three cases sepa- 
rately. 

1. First consider the case when f(m) decays faster than exponential for 
large m. In this case, p is allowed to take any value in [—00, 00] and 
p(p) is monotonic in p, ranging from 00 to 0. Thus, for a given p, 
one can always find a suitable p to satisfy Eq. (P) and there is no 
condensation. 

2. We next consider the case when f(m) decays slower than m~ 2 for large 
m. In this case, the allowed range of p is [0,oo], in order that the 
integrals in © converge. The function p(/i) diverges as p — > 0, and 
then decreases monotonically to p{p) — > as p — > 00. Once again, for 
any given p, there is always a solution of p from Eq. (jSJ and there is 
no condensation. 
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3. Finally consider the case when f(m) decays (a) slower than exponential 
but (b) faster than m~ 2 . In this case, the allowed range of p is [0, oo] 
over which p(p) is a monotonically decreasing function of p. However, 
the crucial difference here from the two previous cases is that p(0) is 
finite. Indeed, p(0) sets the critical density p c = p(0). When the given 
p < p c , a positive p satisfying (0) exists. However, when p > p c , there 
is no real solution to ©, signalling a condensation transition. In this 
phase, the extra mass (p — p c )L forms a condensate. 

A natural choice of f(m) satisfying the criterion of condensation in 3. 
above is 

f(m) ~ Am~ J with 7 > 2. (7) 

For the majority of this paper, we stay with the choice of /(m)'s in (J7J) and 
set, without loss of generality, J °° /(m)dm = 1. For brevity, we also assume 
that 7 is non-integer although our analysis can be easily extended to the 
integer case. In section some other forms for f(m) will be explored. 
The single-site probability distribution p(m) is given by 

p(m) = f(m)e~ fim for p < p c . (8) 

which exhibits a characteristic mass m* = 1/p. As the critical density is 
approached p decreases to zero and the characteristic mass diverges giving 
rise to a power law p(m) ~ f{m) at the critical density. For p > p c the grand 
canonical approach can not be used to determine p{m). 

This last point is in contrast to usual Bose-Einstein condensation where 
one can work in the GCE even in the condensed phase. This is done by 
letting the p tend to zero as 1/V where V is the volume of the system i.e. 
any density of bosons can be achieved by carefully letting p \ in a way 
dependent on system size. However, in the present case letting p — > in (JBJ) 
always results in a fixed critical density p c . 

3 Interpretation as sums of random variables 

Before proceeding with a detailed analysis of the canonical partition function 
((21) it is instructive to discuss the problem from the perspective of sums of 
random variables. 

First note that if f(m) is normalised then Ylf =1 is the probability 

that L independent and identically distributed (iid) positive random vari- 
ables, each drawn from a distribution f(m), take the values mi,m2, . . .mi 
which we denote by m. Moreover, the partition function Z(M, L) in Q is 
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precisely the probability that the sum of the random variables is equal to M. 
Equivalently one can think of an ensemble the configurations of which are 
defined as masses m whose sum is M, with the weights of the configurations 
being Yli=i /( m «)- The question of condensation then reduces to the study 
of the statistical properties of the largest of the L random variables. 
Let us define the moments of f(m) as 

POO 

/ik — dmm k f(m) . (9) 
Jo 

If the mean /zi exists (e.g., (|IJ)) and is such that L/ii > M, we expect the 
ensemble to be dominated by configurations where m« = 0(1) Vi. However, 
if L/i! < M then we expect the ensemble to be dominated by configurations 
where L — 1 random variables are 0(1) and one is O(M). Thus, for f{m) ~ 
m~ 7 and 7 > 2 we expect condensation at a critical density p c = H\. This 
recovers the results given in sectional 

We can also obtain the asymptotic behaviour of the partition function 
(|2*j) by noting that the large random variable could be any of the L possible 
ones and that its probability ~ A(M — M c )~ 7 . Therefore we expect in the 
condensed phase 

Z(M, L) ~ AL/{M - M c y . (10) 

It turns out that this simple argument gives the correct asymptotic behavior 
to be derived in section El (see Eqs. (|HU|) . (jl03j0 . 

Let us now consider the fluctuations of the condensate in the case 7 > 2. 
Since we have the constraint of total mass M we can equally well consider 
the fluctuations of the total fluid mass which is essentially the sum of L — 1 
random variables drawn from f(m): 

L-l 

M f = Y,™*- (11) 

i=l 

Now if 7 > 3, then /i 2 is finite and 

A = y/m-nl (12) 

is the well defined width of / (m). The standard results of the Central Limit 
Theorem applies to this fluid component, with 

AM, = [(M})-(M f f] l/2 (13) 
z> [(/i 2 -/i?)L] 1/2 for 7 >3 (14) 
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By constraint, the fluctuation of the condensate will be controlled by A\/L. 

However if 2 < 7 < 3, then p 2 does not exist. Instead, AM/ ~ Ylii m 1 
will be dominated by the largest of the L — 1 random variables in the fluid. 
Therefore 

AM/ = O(L^) for 2 < 7 < 3 . (15) 

In this case the fluctuations of the condensate are large AM/ > 0(L 1/2 ) and 
we refer to this as an anomalous condensate. 

In fact, using the sums of random variables interpretation we obtain very 
quickly the scaling behaviour of our main results. In the following sections 
we will obtain more precise and detailed results. 



4 Nature of the Condensate: Summary of 
Results 

The GCE analysis of section 121 correctly predicts the criterion for condensa- 
tion and even the critical density p c = J °° mf(m)dm (whenever condensa- 
tion occurs), but provides little insight into the condensed phase itself where 
p > p c . In this section we explore the condensed phase in detail by staying 
within the framework of a 'canonical ensemble' and analyzing the single site 
mass distribution 

p(m) = J dm 2 ....dmx J P(m, m 2 , ■ ■ ■ , m^S rrij + m — M^j , (16) 

in a finite system of size L and total mass M = pL. Note that p(m) de- 
pends on L, though we have suppressed its L dependence just for notational 
simplicity. Using Eqs. (JTJ) and (j2J, we have 

, , „. S Z(M -m,L-l) , . 

« = Z(M,L) • < 17 > 

The rest of this section is devoted to the analysis of p(m) in (|T7jl with /(to) 
given by (|7jl. We have thus two parameters 7 and p. Our goal is to show how 
the condensation manifests itself in different behaviors of p(m) in different 
regions of the (p — 7) plane. We will show that for 7 > 2, there is a critical 
curve p c (l) in the (p — 7) plane that separates a fluid phase (for p < p c (j)) 
from a condensed phase (for p > p c (j))- In the fluid phase the mass dis- 
tribution decays exponentially for large m, p(m) ~ exp[— m/m*} where the 
characteristic mass m* increases with increasing density and diverges as the 
density approaches its critical value p c from below. At p = p c the distribution 
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decays as a power law, p{m) ~ m~ 7 for large m. For p > p c , the distribution, 
in addition to the power law decaying part, develops an additional bump, 
representing the condensate, centred around the "excess" mass: 

M ex = M- Pc L. (18) 

Furthermore, by our analysis within the canonical ensemble, we show that 
even inside the condensed phase (p > p c (l)), there are two types of behaviors 
of the condensate depending on the value of 7. For 2 < 7 < 3, the condensate 
is characterized by anomalous non-gaussian fluctuations whereas for 7 > 3, 
the condensate has gaussian fluctuations. This leads to a rich phase diagram 
in the (p — 7) plane, a schematic picture of which is presented in Fig. (|T|). 

20 n . 1 . 1 



ANOMALOUS 
CONDENSATE 



NORMAL 
CONDENSATE 



FLUID PHASE 



Figure 1: Schematic phase diagram in the p-7 plane. 
To proceed, we take the Laplace transform of (j2J): 

POO 

/ Z(M,L)e- sM dM = [g(s)] L , (19) 
Jo 

where 

/■oo 

g (s) = / f(m)e- sm dm . (20) 
Jo 

The main challenge is to invert (fT9"j) for a given f(m) to compute Z(M,L) 
then exploit its behavior via (fTTj) to analyse the single site mass distribution 
p(m). It turns out that in certain cases, to be discussed in sectional one 
can invert (|T9"j) exactly. However, for general f(m) we rely on an asymptotic 
analysis to be presented in section |3 Before presenting the details of these 
calculations we give a summary of our main results for p(m) 
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4.1 Summary of Results 

Fluid phase p < p c : In this case one finds (see I68|69|) 

p{m) ~ f(m) e~ m/m * for 1 « m « M (21) 

where the characteristic mass m* diverges p approaches p c from below as 
(p - p c )~ l for 7 > 3 and (p - p c )" 1/(7 " 2) for 2 < 7 < 3. 

Condensed Phase p > p c : In this case one finds 

p(rn) ~ f{m) for 1 < m < O(L) (22) 

p(m) ~ /(m) — — for m = xM ra where < x < 1 (23) 

p(m) ~ Pcond(^) for m ~ M ex (24) 

Here p con d is the piece of p{m) which describes the condensate: Centred on 
the excess M ex and with integral being equal to 1/L, it takes on two distinct 
forms (see I90I105|) according to whether 7 < 3 or 7 > 3. For 2 < 7 < 3 



Pcond 



(m) ~ L-^^- l) V 7 



m — M e . 



(25) 



where Vy(z) is given exactly by Eq. (|77jl with asymptotic forms in Eqs. (I7S1 
l80j) . Thus the shape of the condensate bump is non-gaussian for 2 < 7 < 3 
and we refer to this as an 'anomalous' condensate. On the other hand, for 

7 > 3 

Pcond (m) ~ - T =L= e -(m-M^f/2A*L for |m _ ^ <<; 0(L 2/ 3)> (26) 

v 27rA 2 L ,:i 

i-e. Pcond( m ) is gaussian on the scale \m — M ex \ -C 0(L 2//3 ), but, far to the 
left of the peak, p{m) decays as a power law (see llOfij) . 

Critical density p = p c : In this case one finds (see ll()8pi()9jl that 

p(m) oc f(m)Vy (m/L 1 /^) , for 2< 7 <3 (27) 
p(m) oc f(m) e - m2/2A2L 7 >3. (28) 

where the scaling function Vy(z) is given by Ij77j) . Thus at criticality p(m) 
decays as a power law m~ 7 for large m which is cut-off by a finite size scaling 
function and the cut-off mass scales as 

m cut _ ofr ~ ^V(7-i) for 2 < 7 < 3 (29) 
~ L 1/2 for 7 > 3. (30) 
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5 Condensation in the Canonical Ensemble: 
Exactly Solvable Cases 

Before proceeding to derive the results of section 14.11 in the general case, it 
is useful to work out cases where both Z(M, L) and p(m) can be obtained 
in closed form. We consider two examples here. In the first example there 
is a genuine condensation transition, whereas in the second one has only a 
pseudocondensation. 

Example I: Making the choice 

/M = ^=4^ e " 1/m ' ( 31 ) 



for all m, yields the large m behavior (JZJ) with 7 = 5/2 and A = Ij^/H. To 
calculate the Laplace transform g(s) of (|3~Tj) we make use of the following 
identity, 3.471.9 from Gradshtyn and Ryzhik[ 



J" x^e-^-^dx = 2 (^j ' K " ( 2 V^7) > (32) 

where K v (x) is the modified Bessel function. Choosing (3 = 1, j = s, v = 
—3/2 and using the fact that K- V (x) = K v {x) one gets 

/>oo 

/ ,• 5 \ 1 ' sx dx = 2s 3 / 4 K 3/2 (2V5) • (33) 
Jo 

Next, we use the expression 

1T(1 + z 



^w = V^ z e "■ (34) 

Substituting this on the r.h.s of Eq. (|3^j) and simplifying we obtain 

g{s) = -== / m- 5/2 e- 1/m - sm dm = (1 + 2^)e~ 2 ^ s . (35) 
V 71 " Jo 

We then substitute this result on the r.h.s. of Eq. (fT9"j) and expand the r.h.s 
to obtain 



^ Z{M, L)e~ sAI dM = £ ( L \ (4s 
^° fc=0 ^ ' 



i fc/ V 2L ^. (36) 



To proceed we need to invert the Laplace transform in Eq. (|36j) . For this, 
we require the inverse Laplace transform, L~ x ^ fe / 2 e _av/ *] with a = 2L. This 
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is done by again using the identity in Eq. (J32*j) . Using u — — 1/2, (3 — a 2 /4, 
7 = s and 

K 1/2 (z) = ^JTze~ z (37) 

we get 

' X 2 '^ — - V^-a^5 ( 38 ) 



x^ 2 e- a ^ x e- sx dx 
This provides us with the identity, 

c: 1 



a e 



-a 2 /4x 



1 d 



da 



-a 2 /Ax 



(39) 



where x is the argument of the inverse Laplace transform. Next we differen- 
tiate both sides of Eq. (J39|) k times with respect to a. This gives, 

fc+i 



c- 1 



s k/2 e ~a^/s 



(_l)fc-l / Q 

tJtxx \da 

(^)- {k+l)/2 ^a 2 ^ 



-a 2 /Ax 



Hk+i 



(40) 



where Hk(x) is the Hermite polynomial of degree k and argument x and we 
have used the Rodrigues' representation 



Hdx) 



(41) 



Using the result (jjOj) , we then invert the Laplace transform in Eq. ()3fij) to 
obtain 



Z(M,L) 



1 e~ L2 l M 
2^ M 



fc=0 



(42) 



We still need to perform the sum in Eq. (J42j) . To do this, we use the 
following identity which is proved in the appendix-A, 



£ (f) b~ k H k (x) = b- L H L (x + b/2), 

k=0 ^ ' 



(43) 



where b is any constant. Now, consider the sum 



E 

fc=0 



L\ H k+1 (x) 
k b k 



E 

fc=0 



-k 



2x H k (x) 



dH k (x) 
dx 



d 



2x b' L H L {x + 6/2) - — [b~ L H L {x + b/2)] 



2x b~ L H L (x + b/2) - 2L b~ L H L ^(x + b/2). (44) 



-L 
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In going from 1st to the second line of Eq. (|44|) . we have used the well 
known identity, H k+ i(x) = 2xH k (x) — dH k /dx, in going from 2nd to the 
3rd line we have used the identity in Eq. ()43)1 and in going from 3rd to the 
4th line we have used another well known identity dH k (x)/dx = 2kH k -i(x). 
Substituting this result on the r.h.s. of Eq. (|42jl and simplifying we obtain 
our final formula, valid for all M and L, 



Z(M,L) = J- M -^ 2 e- L2 / M 

'IT 



V l^fM ) \ 2^M 

(45) 

As a check for consistency, one can easily verify that for L — 1, using H (x) = 
1 and Hi{x) = 2x, the formula in Eq. (j43j) reduces to 

= 4 jj^S«- ,/ " = /(*). (46) 



as it should be for L = 1 from the definition in Eq. (J2J). 

Since f(m) ~ m~ 5//2 for large m, we expect from the GCE analysis in 
section-3 that in this case there will be a condensation transition at the 
critical density, p c = /?(// = 0), where p{p) is given by Eq. (jHJ). Now, 
p(/i = 0) can be obtained from the Laplace transform: p(p = 0) = —g'(0), 
where g'(s) = dg/ds. Using Eq. (j33|) we get 

Pc = 2 . (47) 

To see how the condensation transition manifests itself in the single site mass 
distribution function p(m), we calculate p(m) explicitly by substituting in Eq. 
()17|) the expression for f(m) from Eq. (|3*T|) and that of Z(M,L) from Eq. 
(j45|) . Using Mathematica, we plot this explicit expression of p(m) for L = 100 
for three densities, p = 1 (subcritical), p = p c = 2 (critical) and p = 6 
(supercritical) in Fig. |21 The transition from the fluid phase (p < p c = 2) to 
the condensed phase (p > p c = 2) is clearly seen: the condensate appears as 
an additional bump in p(m) near its tail for the supercritical case p = 6. 

The exact solution detailed above is instructive. It confirms the grand 
canonical prediction that the signature of the condensation transition is man- 
ifest in the behaviour of the single site mass distribution p(m) near its tail, as 
evident in Fig. (J2J). The asymptotic decay of p(m) for large m is different for 
p < p c , p = p c and p > p c . For p < p c , p(m) decays exponentially for large 
m. For the critical case p = p c , p(m) has a power law tail: p(m) ~ m~ 5 ^ 2 
for large m. As p increases beyond p c , the critical power law part of p(m) 
does not change, but the extra mass (p — p c )L accumulates on the condensate 
which shows up as an additional bump beyond the power law tail of p(m). 
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Figure 2: The distribution p(m) vs. m for the exactly solvable case, plotted 
using Mathematica for L = 100 and p = 1 (subcritical), p = p c = 2 (critical) 
and p = 6 (supercritical). The condensate shows up as an additional bump 
near the tail of p(m) in the supercritical case. 

Thus for p > p c , p(m) has two parts: a power law decaying part followed 
by a bump, indicating that the condensate coexists with a background that 
behaves as a critical fluid. Even though this exact solution was obtained 
only for a specific /(to) (pTT|). it confirms the generic picture of section |4~T1 
for the behaviour in the tail of p(m) at subcritical, critical and the super- 
critical densities respectively. Our goal in the next section is to analyse the 
large m behaviour of p(m) for these three cases for a general f(m) allowing 
a condensation transition, as in Eq. (JZJ). 

Example II: It is also instructive to solve exactly where f(m) 

decays more slowly than 1/m 2 and condensation does not occur. This leads 
us to our second exactly solvable example where we find that even though 
there is no genuine condensation, a pseudocondensate appears nevertheless. 

We consider 





for which, using (jHZJ), one finds 




(49) 



Inverting (JT!5j), using (J39"j) . yields 




(50) 
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Figure 3: The scaling function Y(x) vs. x in Eq. (|52|) . for (p — 1, (p — (p c — 
3.88562 and p = 6. The pseudocondensate shows up as a bump near the tail 
of Y(x) for (p > p c . 



and from ()17|) we find 

1 L — 1 / M \ 3 / 2 
p(m) = — — exp(-(L-l) 2 j[M -m) + L 2 /M -l/m) 



7T 1 / 2 L \m(M — m 

(51) 

A natural scaling regime emerges when M = ipL 2 ,m = xM. Note that 
in this regime the density is superextensive. Taking the limit of L large we 
obtain p(x M) = 7r~ 1/2 M~ 3/2 Y(x) where 

Y(x) = l x)Y /i exp[-x/(y(l - x))\ . (52) 

A bump in the tail of this distribution emerges when there are two turning 
points of Y(x) at positive x. This occurs when ip > ip c where ip c = 2 + 
Ay/2/3 = 3.88562 . . . (see Fig.EJ). As we shall discuss later in section IT2l this 
bump is not a true condensate but rather corresponds to what we shall refer 
to 'pseudocondensate' . 
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6 Condensation in the Canonical Ensemble: 
General Case 

We now proceed to analyze the case of a general f(m) characterized by a 
power law tail with exponent 7 > 2 as in Eq. (J7J). As already noted the 
Laplace transform 

POO 

g(s) = / f(m)e- sm dm (53) 
Jo 

plays a crucial role in our analysis. We start by formally inverting the Laplace 
transform in Eq. (j!9|) using the Bromwich formula, 

pso+ico j 

Z(M,L)= — exp [L (In g(s)+ps)] (54) 

Jso-ioo ^ m 

where we have used M = pL. Similarly, one has 

r-so+ico ^ 

Z(M — m,L) = / — : exp [L (\ng(s) + ps) — sm] . (55) 
Jso—ioc 27ri 

In Eqs. (I54j) and 1)55)1 . the contour in the complex s plane parallels the 
imaginary axis with its real part, Sq, to the right of all singularities of the 
integrand. Since f(m < 0) = 0, the integrand is analytic in the right half 
plane. Therefore, so can assume any non-negative value. Meanwhile, for / 
given by (|ZJ), s = is a branch point singularity. As we shall see, in the 
subcritical case p < p c there exists a saddle point at positive s and so can 
be chosen to be this saddle point (see Fig. (j3|). One can then evaluate the 
Bromwich integral in Eq. by the saddle point method for large L. As 
the density p approaches its critical value p c from below, the saddle point 
so moves towards 0. Thus, in the critical (p = p c ) and in the supercritical 
(p > p c ), one can no longer evaluate the Bromwich integral by the saddle 
point method and an alternative approach must be used. 

Let us first evaluate the integral in Eq. ()54|) in the limit L — ► 00 by the 
saddle point method, assuming it exists. It is convenient to define 

h(s) = ps + In g(s). (56) 

Then the saddle point equation, h'(so) = 0, is 

P = -g'(so)/g(s ) } (57) 

which is precisely equivalent to GCE approach in Eq. (jSJ on identifying 
So = /i. Thus the saddle point evaluation works provided there exists a 
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Figure 4: The contour parallel to the imaginary axis passing through the 
saddle point s in the complex s plane. 

solution s > of Eq. (|57j). As clarified in the GCE analysis, a solution 
exists when p < p c where 

Pc = -9'(0)/g(0) = -g'(0). (58) 

We have assumed that f(m) is normalized to unity so that g(0) = 1. Also, 
when the saddle point So > exists, one simply has 

Similarly when m <C O(L), from Eq. (j^} one gets Z(M—m, L) ~ Z(M, L)e~ s 
Substituting in Eq. (fTTjl we get, for p < p c and l<Cm-C O(L), 

p(m) ~ /(m)e- Som , (60) 

thus recovering the GCE result upon identifying so = This shows that 
for p < p c , the canonical and the grand canonical approach are equivalent 
and the mass distribution p(m) decays exponentially for large m with a 
characteristic mass m* = 1/sq. 

As p approaches the critical value p c from below, so approaches and 
hence the characteristic mass m* diverges. To see how the asymptotic be- 
haviour of p(m) changes as one increases p through its critical value, we will 
now focus only in the vicinity of the critical point. Near the critical point, i.e, 
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when \p — p c \ is small, the most important contribution to the integral in Eq. 
(|54|) comes from the small s region. Hence, to obtain the leading behaviour 
for large L it is sufficient to consider only the small s behavior of g(s) in Eq. 
(jHIJ). For f(m) in Eq. (J2J) with a noninteger 7 > 2, one can expand, quite 
generally, the Laplace transform g(s) in Eq. (|53j) for small s, as 



n-1 

^) = E(- 1 ) fc t sfc+fo7_1 + --- (61) 



fc=0 

Here n = int[7], pk is the k th moment of f(m) (which exists for k < n). We 
assume that f(m) is normalized to unity, i.e., po = 1. Note also that from 
Eq. dSHJ) it follows that 

Pi = Pc- (62) 

The term 6s 7_1 in Eq. (jfilj) is the leading singular term. For /(m) in Eq. (jJJ) 
with a noninteger 7, one can show (see appendix-B) that 

6 = Ar(l- 7 ), (63) 

a relation that we will use later. Note that the expansion in Eq. (JBTJl is valid 
only for noninteger 7. For integer 7, one can carry out a similar analysis and 
it is easy to show that the leading singular term has a logarithmic correction: 

7-2 

g (s) = Y^(-l) k ^s k + cs^ 1 \ns + ... for 7 = 2,3,... (64) 

fc=0 

where 



k\ 

(-1)^4 



(65) 



(7-1)! 

Substituting the small s expansion of g(s) from Eq. (jBTJl in Eq. (jB^j) and 
keeping only the two leading order terms explicitly one gets 

h(s) = (p- p c )s + 6s 7_1 + ... for2<7<3 

= (p-p c )s- As 2 log(s)/2 + ... for 7 = 3 

= (p-p c )s + A 2 s 2 /2 + --- + fe 7 - 1 + ... for 7 >3 (66) 

where A is defined in Eq. f(12jh The special role of 7 = 3 is now clear. The 
next-to-leading term in h(s) is s 7_1 for 2 < 7 < 3 and s 2 for 7 > 3. 

Below we substitute the small s expansion of h(s) from Eq. (|66j) in the 
integral in Eq. (p>4"|) and analyse its leading asymptotic behaviour for large L 
for the three cases p < p c , p > p c and p = p c separately. Since we are using 
the small s expansion, this analysis is valid only in the vicinity of the critical 
point, i.e., when \p — p c \ is small. 
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6.1 Subcritical Case: p < p c 



For p < p c , it follows from Eq. (jfifijl that h(s) has a saddle at so > 0. This is 
equivalent to saying that Eq. (|57jl has a solution which, to leading order in 
(pc - p), is given by 

so * for2< 7 <3 
.6(7 - !). 

~ -{p c - p) / A\og{p c - p) for 7 = 3 

~ (p c - p)/A 2 for 7 > 3. (67) 

Consequently, it follows from Eq. (fTT]). as well as from Eq. (jSOjl . that the 
single site mass distribution p(m) behaves for l<m< 0(L) as 

v ' nV v ' 

where sq is given by Eq. (|fi7j) for small (p c — p). Thus the characteristic mass 
m* = 1/so diverges as a power law as p approaches p c from below, 

m* ~ (p c - p )- 1 /(7~2) for 2 < 7 < 3 

~ - log(p c - p)l (p c ~ p) for 7 = 3 

~ (p c - pY 1 for 7 > 3. (69) 

The partition function in Eq. (|59jh which we will require later, may also 
be calculated and behaves asymptotically as 

Z(M = pL,L) ~ exp [-LB{p c - p)^- t]/ ^- 2) ] for 2 < 7 < 3 
~ exp [L(p c - p) 2 /2A\og(p c - p)] for 7 = 3 
~ exp [-L(p c - p) 2 /2A 2 ] for 7 > 3, (70) 

where B = (7 - 2) 6~V(7-2) ^ _ 1 ^_( 7 _i)/( 7 _2)_ 



6.2 Supercritical case: p > p c 

For the supercritical regime (p > p c ), there is no solution to (jf>7j) on the 
positive real axis and more care is needed to find the asymptotic form of 
Z(M, L). One can see this clearly in the vicinity of the critical point where 
one can use the small s expansion of h(s) in Eq. evidently for p > p c , 

h(s) does not have a minimum at any positive Sq. In the L — > 00 limit and 
with (p— p c ) small, one can, however, develop a scaling analysis by identifying 
the different scaling regimes and calculating the corresponding scaling forms 
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for Z(M, L) and hence that of p(m). The main result of this subsection is to 
obtain an exact asymptotic formula describing the shape of the condensate 
bump. 

Substituting Eq. fDJ in Eqs. and (jHKjl . we get 

Z(M,L) ~ ^l(Pc-p) (71) 
Z(M-m,L) ~ ^((m-ikQ/L), (72) 

where M ex = (p — p c )L is the excess mass and 

W L (y) = J l °° ^exp [L(-ys + C(s))] . (73) 

In Eq. (J75j) the function C(s) represents the nonanalytic correction to the 
linear term in h(s) in Eq. ([66)1. i.e., 

C(s) = 6s 7 " 1 for 2 < 7 < 3 

= -As 2 log(s)/2 for 7 = 3 

= AV/2 + ■ ■ • + 6s 7 " 1 for 7 > 3 (74) 

By including the first nonanalytic term in C(s) we will obtain the leading 
large L asymptotic behaviour. 

Knowing the function Wl{v) and using Eqs. (J7T|) and (J72j) . one can 
rewrite the mass distribution in Eq. ()17|) as 

P(m) = /W ^«7 M 'f £) . (75) 

(Pc - P) 

All crucial information about the condensate 'bump' around m = M ex = 
(p — p c )L is thus encoded in the asymptotic behavior of Wl(v) defined in Eq. 
(|75|). Below we analyze the asymptotics of Wl(v) and hence that of p(m) in 
the two cases 2 < 7 < 3 and 7 > 3 separately. 

Case-I 2 < 7 < 3 : In this case, we substitute C(s) = 6s 7_1 in Eq. 
(|73|) and rescale s — > L -1 '^ -1 ^ to rewrite Eq. (|73|) in the scaling form 

W^(y) = L- 1 ^- 1 ^ [L^ a )/(r-i) y ] , (76) 

where the scaling function is 

^L^'* (77) 
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We were not able to perform the complex integral in Eq. (|77jl in closed form. 
However, its asymptotic behaviors can be worked out, the details of which 
are relegated to appendix-C. We find, 

Vy{z) ~ A\z\^ asz^-oo (78) 
= c at z = (79) 

~ Cl2 (3-7)/2(7-2) e - C2 ^-D/(^) aaz ^ QQ (80) 

where A is the amplitude in Eq. (|7J), 5 = A r(l — 7) (as in Eq. (|63j)) and the 

constants c , ci and c 2 are given as 

c = 6-i/(^i)/[( 7 _i)r(( 7 -2)/( 7 -l))] (81) 
Cl = [27r( 7 -2)(6(7-l)) 1/( ^ 2) ] _1/2 (82) 
c 2 = (7-2)/(7-l)(K7-l)) 1/(7 ^ 2) . (83) 

Evidently the scaling function V 1 (z) is highly asymmetric, has an al- 
gebraic decay |z| -7 as z — > —00 and decays extremely fast (faster than a 
gaussian) as z — > 00. To plot the full function V^(z) in Eq. f77j). it is useful 
to first transform the integration in the complex s plane in Eq. (J77|) to a real 
integral. This can be achieved by making a change of variable s = e ±m ' 2 y 
respectively for the upper and the lower half of the imaginary axis in the 
integral in Eq. (J77j) and then simplifying. One obtains 

1 f°° 

V y (z) = - dye' 03 ^ cos [6cos(vr7/2)?/ 7 - 1 + yz] , (84) 
71 Jo 

where C3 = — 6sin(7T7/2) > for 2 < 7 < 3. The real integral in Eq. (jSlj) can 
be easily evaluated numerically. A plot of this function for 7 = 5/2, where 
the integral in Eq. (jSljl was performed numerically using Mathematica, is 
shown in Fig. (JHJ). The dashed lines at the two tails show the agreement 
with the asymptotic analytical forms in Eqs. (fT8*|) and (|8()j) . 

Armed with knowledge of the asymptotic behaviors of Wl{v) we are now 
ready to calculate the mass distribution p(m) from Eq. J75J). Let us first 
evaluate the denominator in Eq. (|75|) using the scaling form of Wi(y) in Eq. 

(Hi, 

Z(M, L) ~ W L ( Pc -p)= L _1/(7_1 V 7 [-L^- 2 )/^- 1 ) (p - p c )] . (85) 

Since p > p c , the argument of the scaling function in Eq. (}85|) is a large 
negative number as L — > 00, provided p — p c = M ex /L > is kept fixed. One 
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Figure 5: Plot of the scaling function Vy(z) for 7 = 5/2 obtained by nu- 
merically integrating the integral in Eq. (j84|) . The dashed lines show the 
asymptotic tails obtained analytically in Eqs. (J7%|) and (JHUj) . 



can then use the asymptotic tail of Vy(z) as given in Eq. (|7Hj) and we get, in 
the regime where M ex ~ O(L), 



Z(M,L)~W L (p c -p) 



A 



AL 



L^' 1 (p - PcP M ex y 
Collecting all these results together in Eq. (|75j) one then finds 



p(m 



)~ZM( p _ pc )7 L 7(7-2)/(7-Dy 7 



A 



m - M e: 
LV7-1 



(86) 



(87) 



Note that the result in Eq. (J%7|) is valid for all m, provided the total mass 
M and and the system size L are both large with Mj L = p fixed. Now one 
can identify various limits of (|87jl . Using (fTSj) one finds 



p(m) ~ f( m ) f° r l<m< O(L) 
p(m) ~ f(m)-, — for m = xM e 

V ^ V 7 (1-X)7 



where < x < 1 (89) 



Thus for m O(L), p(m) is a pure power law but when m becomes extensive 
p(m) begins to deviate from f(m). 

Focusing now only near the condensate where m ~ M ex , one obtains using 
f(m) ~ Am' 1 the asymptotic form of p(m) ~ Pcond(^) near the condensate 
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bump, 

/ v r _ 7 /( 7 _i )T/ \m-M P .. 



LV7-1 



(90) 



where V^(z) is given exactly by Eq. (|77jl with asymptotic forms in Eqs. (JTSJ- 
IHU|) and a full shape as shown in Fig. (|5|). Thus the shape of the condensate 
bump, as shown in Fig. (jSJ) is highly asymmetric and non-gaussian for 2 < 
7 < 3 and we refer to this as the 'anomalous' condensate. Also note that, 
setting £ = (m - M^L' 1 ^- 1 , 

/oo 
df^ 7 (0 (91) 
-oo 

Thus the area under p con d(^) is 0(1/ L) which corresponds to precisely one 
condensate in the system of size L. 

Case-II 7 > 3 : In this case, using the appropriate form for C(s) from 
Eq. (JUJ) in Eq. {Efl we get 

r°° ds 

W L (y) = — exp L(-y S + AV/2 + • ■ • + bs^ 1 ] . (92) 

We will see that Wl(v) has different behaviors for y > and y < 0. Let us 
first consider the case y > 0, where one can again perform a saddle point 
analysis as in the previous subsection and one gets, to leading order in L and 
for yy/L ~ 0(1) 

W L (y) * J-r^ z-^. (93) 
V 2 7T A z L 

Thus the function W^y) has a gaussian decay for positive but small y. A 
careful analysis shows that the gaussian form with the scaling variable yy/L 
in Eq. (J93|) remains valid not just over the range yy/L ~ 0(1), but over a 
wider range: at least up to y ~ 0(L -1 / 3 ). For large y > 0, one can still do a 
saddle point analysis 'a la Eq. (jHUjl and Wi(y) will have a non-gaussian tail 
for large, positive y. 

We next turn to y = — \y\ < 0. In this case, there is no saddle point and 
we have to do the integration along the imaginary axis passing through the 
origin. Again, let us first consider the regime where \y\yL ~ 0(1). In this 
case, we rescale s — > s/y/L in Eq. (|9*2~j) which then becomes 



1 f iaa d s r 
W L (y) = -= — exp -(yy/L)s + AV/2 + • • • + bL& 
y/L J-ioo 27rz L 



•7)/2 s 7-l 



(94) 
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Keeping \y\yL = z fixed and taking the L — ► oo limit, one can drop all the 
subleading terms and the resulting integral is a gaussian one which can be 
simply performed to give for \y\ ~ 0(l/vL), 



W L (y) ~ / -_— e^ 2L / 2A2 , (95) 

the same result as for y ~ 0(L~ l l 2 ) > in Eq. (|93|) . Thus, within the range 
|y| ~ 0(L -1 / 3 ), the function H^d/) is symmetric and gaussian. On the other 
hand, far to the left of the origin where \y\ 3> 0(L -1 / 3 ), this analysis breaks 
down. Let us illustrate the case \y\ ~ 0(1). In that regime of we again 
start from Eq. ()92)) but this time rescale s — > s/L which gives 





MOO 


ds 


L t 


' — ioo 


2tH 




rico 


ds 


~L „ 


' — ioo 


2tH 



W M = Tl ^-^P[\y\s + ^ 2 s 2 /2L + --- + bL 2 -'s^ 1 ] (96) 

(97) 



A 2 s 2 , s 7 " 1 

1 + ^ + --- +6 z^ 



where in going from first to the second line, we have expanded the exponential 
for large L, keeping \y\ fixed. Now, the integral in Eq. ()97j1 can be done 
term by term. The first term gives a delta function S(\y\) which is since 
\y\ > 0. Similarly the second term, which is the second derivative of the 
delta function with respect to \y\, is also zero. In fact, all the analytic terms 
containing integer powers of s will be similarly equal to 0, except the last 
singular term which has s 7 " 1 . Thus one gets for \y\ ~ 0(1) with L — ► oo, 

h f' 00 ds h r°° ds 

W M = L 2S ^ S ' ^ = WL^ L SS ' S7 "' (98) 

The latter integral is done in appendix-C (see Eq. ()138j) ). Using this result 
in Eq. ([9*H)l one obtains 

W L (y) ~ A L7 _, for \y\~0{\) (99) 

Although here we consider \y\ ~ 0(1) it is easy to show the result also holds 
for \y\ ~ 0(L~ a ) where a < 1/2. 

Thus, the large L behavior of the function Wiiy) for 7 > 3 can be 
summarized as follows, 

W L (y) ~ 1 e ~ y2L/2A2 for \y\ < 0(L- 1/3 ) (100) 

v 2 7r A 2 L 

for -j/= ~ 0(1) (101) 
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The function Wl(v) is thus a gaussian near its peak at y — and then as 
non-gaussian tails far away from the peak. On the negative side, this tail is 
algebraic and decays as |y|~ 7 . 

We now compute the asymptotic behavior of the mass distribution p(m) 
in Eq. (fTHJ) for 7 > 3. First, we calculate the partition function in Eq. (fTTf . 
Substituting the results in Eqs. (jl00|) and (jlOlj) in Eq. (|7T|) we find 



Z{M,L) ~ W L (-(p- Pc )) 

e -M ex 2 /2A 2 L 



ioi\M ex \~0{L 2 l*) (102) 



V2 7T A 2 L 

for M e:c ~ O(L). (103) 



Comparing Eqs. ()103j) and (JHHJl one sees that for M ea; > (i.e., M > p c L) 
and in the regime where it is O(L), the partition function Z(M,L) has the 
same asymptotic form for all 7 > 2. Using the results in Eqs. (jlOOj) and 
(|101|) one can similarly evaluate the numerator in Eq. (|75j) . Collecting these 
results, we find that 

t \ t< ^ M - 7 TI/ (m-M ex \ 
p(m) ~ f(m) -jj- W L I J . (104) 

The the asymptotic form of p(m) = p CO nd( m ) near the condensate bump has 
a gaussian peak centered at m — M ex = (p — p c ) L with its width scaling as 
0(\/L), but has a non-gaussian tail far away from the peak. More precisely, 
one gets the following behavior near the peak: 

Pcond (m) ~ - 1 e ~(m-M ex f/2A 2 L for (m _ ^ 0(jL 2/ 3) _ (1Q5) 

On the other hand, for m far less than the peak value of M ex , we have 

p cond (m)~ f(m)(l-m/M ex )-^ for M ex - m ~ O(L). (106) 

The integral of p con( j gives 1/L + 0(L 1-7 ) where the dominant contribution is 
from the gaussian peak. Again this implies a single condensate in the system. 

6.3 Critical Case: p = p c 

In the critical case p — p c , it follows from Eq. (fTHj) 

P W = /H^. dor) 
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We now consider the two cases 2 < 7 < 3 and 7 > 3 separately. 
Case-I 2 < 7 < 3: In this case, using Eq. (f7fi|) we get 

P(m) ^ firn ^^l 1 ^ = ^ V 7 (m/L^) , (108) 

V 7 [\J) Co 

where the constant Cq is given in Eq. (|81|) and the asymptotic behavior 
of the function V y (z) can be read off Eqs. (JTSJ) and dHOJ). The Eq. (fTUHJ) 
thus describes the finite size scaling function associated with the distribution 
p(m): it decays as a power law p(m) ~ f{m) ~ AmT 1 for large m and 
then is cut off by the finite size of the sample and the cut-off mass scales as 
mcut-off-i 1 ^ 1 '. 

Case-II 7 > 3: In this case, using Eq. in Eq. (|107|) we get 

p(m) ~ /(m) e" m2 / 2A2i for m ~ 0(L 2/3 ). (109) 

Thus, in this case, the cut-off function is Gaussian up to m ~ 0(L 2//3 ), even 
though the cut off mass scales as m cut - ff ~ L 1 ^ 2 for all 7 > 3. 



7 Other cases 

So far we have only considered f(m) with asymptotic form ((7j). In this section 
we extend our results to two other choices of f(m). 

7.1 Stretched exponential case 

First we briefly discuss the case 

f(m) ~ Aexp(— cm a ) where a < 1 . (HO) 

This stretched exponential decay clearly fulfils the criterion for condensation 
of section |21 Noting that all moments fi^ of / exist, one finds that the 
subcritical partition function in the vicinity of the critical point behaves 
asymptotically as 

Z{M = pL, L) ~ exp [-L( Pc - p) 2 /2A 2 ] (111) 
and the subcritical p(m) behaves as 

p(m) ~ Aexp(— cm a — Som), (H2) 
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where Sq = (p c — p)/A 2 . Thus in the subcritical regime, there are two 'mass' 
scales. The first one is a natural scale of 0(1) that characterises the stretched 
exponential decay of f(m) itself and it remains of 0(1) even at the critical 
point. The second mass scale 1/sq ~ (p c — p) -1 , that emerges out of collective 
behavior, however diverges as the critical point is approached from below. 
Thus, unlike the power law f(m) discussed in the previous section where 
at the critical point p(m) becomes scale free, for the stretched exponential 
case only one out of the two scales diverges at the critical point. In this 
sense, the condensation transition associated with the stretched exponential 
case is somewhat different from the usual scenario of a second order phase 
transition. 

In the condensed phase p > p c , one can make an analysis qualitatively 
similar to the power-law case with 7 > 3, though the details are somewhat 
different. Omitting details, we find that for p > p c the condensate bump is 
gaussian near its peak as expected 

p cond (m) ~ - 1 e -(m-M ex f/2A*L for (m _ Mgx) „ (m) 

V 2nA 2 L- i 

On the other hand, far to the left of the peak where M ex — m ~ 0(L) we 
have 

Pcond(m) ~ /(m) exp(c [M: x -\m- M ex \ a )) (114) 

7.2 Case 1 < 7 < 2 

As noted in sectional in the case where 

f(m) ~ Am" 7 1 < 7 < 2 (115) 

we do not expect true condensation. However, as we shall now analyse, an in- 
teresting phenomenon may occur at suitably chosen superextensive densities 
where a 'pseudocondensate' bump emerges. 

For ()115j) the corresponding small s expansion of g(s) reads 

g(s) = po + fe 7 " 1 + • • ■ (116) 

where b = AT(1 — 7) and we take as usual po — 1- Now since b < for 
1 < 7 < 2 we always find a solution sq > to the saddle point equation (|H7jl 
which when small is 

r -6(7-l)L ]^ 
so^ ^ (H7) 
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and the saddle point expression for Z(M, L) reads 



Z(M,L) ~ (2tt( 7 - 2) (7 - l)bLs^ 3 ) 1/2 exp 



2-7 f -b( 7 -l)L 

7 - 1 V M^- 1 



) 



1/(2-7) 



(118) 



Let us consider the regime where 



M = tpL 



1/(7-1) 



in 



xM 



(119) 



where (p and 2 are fixed as L — ► 00. Using (JT7j) we find 



p(x) ~ T^-X- 7 (l - x )-(3-7)/2(2-7) exp \_ C ((1 - x)- (7 " 1)/(2 - 7) - 1)] (120) 



As ip is increased a bump in the tail of p(x) emerges corresponding to two real 
turning points of p(x). However we argue that this does not correspond to a 
true condensate for the following reasons. Firstly, as M increases, Z(M, L) 
will follow the same expression ()118|) thus there can be no true phase tran- 
sition since Z(M, L) is analytic. Secondly, there is no diverging length scale 
in p{m). Thirdly, the bump in p(x) is broad i.e. extends over a finite range 
of x, unlike a true condensate bump which would be narrow. Therefore we 
call this bump a 'pseudocondensate' as it is really an extension of the fluid 
phase, rather than a condensate coexisting with the fluid. 

To understand why ()119j) is the natural scale for M, let us consider the 
sum of L random variables each drawn from a distribution ()115j) . We expect 
the largest of the L random variables drawn from f(m) to be ©(L 1 ^ 7-1 )) 
and the total mass to also be 0(L l l^~ l >). Therefore when M is on the scale 
f)l 19j) the constraint of fixed total mass leads to non-trivial p(m). 

8 Discussion 

In this work we have presented an analysis of the partition function (|2*|) of a 
factorised steady state within the canonical ensemble. The analysis has re- 
vealed the structure of the condensate as summarised in section E~T1 The two 
types of condensate occurring when f(m) ~ vrC 1 with 2<7<3or7>3 
correspond to condensates with anomalous (non-gaussian) and gaussian fluc- 
tuations respectively and each has a distinct scaling function for the shape 
of the condensate. In both cases the nonequilibrium phase transition to the 



where 




(121) 
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condensed phase is continuous in the sense that the characteristic mass of 
the exponential site mass distribution of the fluid phase diverges as the tran- 
sition is approached from below. We have also analysed the case of stretched 
exponential f(m) where condensation occurs but the transition is somewhat 
different: even though there is a diverging scale one gets a stretched expo- 
nential rather than power law distribution p(m) at the critical point and in 
the fluid component of the condensed phase. Finally we have shown that in 
the case where f(m) ~ m -7 where 1 < 7 < 2 a pseudocondensate appears 
at superextensive critical density, but there is no phase transition. 

Our results may easily be generalised to the case of discrete mass, ex- 
emplified by the ZRP where the occupation of each site is 0,1,2... and the 
total mass M is a positive integer. In this case the relevant expression for 
the canonical partition function is 

Z(M) L) = / 2^ Z ~ iM+1) [F{Z)]L ' (122) 

where the integral is around a closed contour about the origin in the complex 
z plane and 

00 

F(z) = £ * m f( m ) ■ (!23) 

The phase transition occurs when the saddle point of the z integral reaches 
the radius of convergence of (|12Hjl . For f(m) of the form (|ZJ) the radius of 
convergence is z — 1, therefore we let z = (1 — u) and expand F(z) for u 
small as 

rzi (- v )kp(k) 

k=0 

where r is the integer part of 7 and F^ = 4-% F(z)\ z=1 The second term in 
()124j) is the leading singular part and the coefficient b is again b = AT(1 — 7). 
The general analysis and results then follow analogously to that of section El 
The structure of the condensate also has implications for the dynamics 
within the steady state. In the steady state a site must be randomly selected 
to hold the condensate thus the translational symmetry is broken. However 
on any finite lattice there will be a timescale tl over which the condensate 
dissolves and reforms on another spontaneously selected site. In systems with 
symmetry breaking, the 'flip time' tl [H] is of interest. A recent study [12] 
has shown that for the ZRP with asymmetric mass transfer and dynamics 
which yield f(m) ~ m -7 the flip time grows as ~ L 7 . This result was 
found numerically and also from a simple effective desciption for the conden- 
sate dynamics. That effective description requires a structure for p(m), in 
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particular a dip to the left of the condensate bump, which is borne out by 
the exact results presented here. 

Finally we mention that it would be interesting to analyse condensation 
in steady states that are not factorised, such as in the chipping model of [21] . 
Generally, such steady states exhibit stronger correlations than a factorised 
one. In such cases a factorised steady state is often used as a 'mean-field' ap- 
proximation in the sense that some of the correlations in the true steady state 
are ignored. Condensation is known to occur in some cases of non-factorised 
steady states [21]) on the other hand there are some examples where the mean 
field approximation predicts a condensation transition but in the true steady 
state there is none j2H] ■ It is of importance to further investigate these issues. 
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A Proof of an addition theorem for Hermite 
Polynomials (1431) 

In this appendix, we prove the following identity, 

£ (f) b~ k H k (x) = b- L H L (x + b/2), (125) 

fc=0 ^ ' 

The generating function of the Hermite polynomials is well known, 

Y,u H ^ = e ~ t2+2tx - ( 126 ) 



k=0 



One also has trivially 



X1 t m b 



E^ = «* (127) 



ml 

m=0 



Multiplying Eqs. (fT2fi|) and fTTFj) one gets 



00 -tm+kurn 00 j.n 

£ -lJr E M = e- ,2+2 ' (l+V2 » = E + «/2). (128) 



m\k\ ^— ' n! 

k,m=0 n=0 



Matching powers of t L on both sides of Eq. ()128j) one arrives at the result in 
Eq. flZfl ■ 
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B Proof of expansion ( I61U63I) 

In this appendix we show that for f(m) in Eq. (jZJ with a noninteger 7 > 2, 
the coefficient 6 in the singular term in the small s expansion of the Laplace 
transform g(s) in Eq. ()61|) is related simply to the amplitude A of the power 
law tail of f(m) via 

6 = Ar(l- 7 ). (129) 

Note that n in Eq. ()61|) is simply n = int[7]. Differentiating Eq. (|61)1 n 
times with respect to s and using the definition of g(s) in Eq. (|H3|) . we get 

(-l) n / e- sm m n /(m)dm = 5-^-^ 1 + .... (130) 
Jo r (7 - n) 

Let us denote the l.h.s. of Eq. ()130|) by I(s). Making a change of variable 
y = sm in the integral on the l.h.s of Eq. ()130|) . yields 

POO 

I(s) = (-I)"*""- 1 / e-yy n f(y/s)dy. (131) 
./o 

Next we take the s — > limit in Eq. 1)131)1 . In that limit, the leading 
contribution to the integral will come from the region where the argument 
y/s of f(m) is large, so that one can use f(y/s) ~ Ai^s/y) 1 . Substituting 
this expression and performing the resulting integral, gives, to leading order 
in s, 

I(s) = (-l)MT(n + 1 - 7)s 7 " n - 1 + . . . (132) 
Comparing the leading terms on the l.h.s. and the r.h.s. of Eq. (1130)1 we get 

(-D n 

r(n+l-7)r( T -n)A (133) 



r(7) 

Simplifying Eq. ()133)) using the well known identity, 

T(x)r(l - x) = vr/sin(7rx) , (134) 
gives the result in Eq. (J129jl . 

C Asymptotic Expansion of Vy(z) 

In this appendix we analyse the scaling function Vy(z) 

VJz) = / e -zs+bs->- where b = _ 7 ) ^ ^ 35 j 

J-ioo 2m 
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near its tails z — > ±00 for all 2 < 7 < 3 and also evaluate the function at 
z = 0. 

The tail z — > —00: We write z = —\z\ in Eq. (|77|) and rescale s — ► 
to write 



1 

N 
1 

Ul 



400 ds_ 

27T? 

-200 1 

iDO ds 



, s +b(s/\z\r~ i 



l + a- 



f'- 1 
vIt — 1 



+ 



(136) 



In going from the first to the second line in Eq. (jl36j) we have expanded 
exp[6(s/|z|) 7_1 ] in a power series. We now perform the integration in Eq. 
(j!36|) term by term. The first term is simply zero for any nonzero \z\. The 
second term is of order 0(|z|™ 7 ). In general, the (n + l)-th term will scale 
as 0(\z\ n ^' y ~ 1 ' +1 ). Thus, for large \z\, the leading asymptotic behavior is 
captured by the second term and one gets 



\z\i 



ds 
2tH 



(137) 



This integral may be performed by wrapping the contour around the branch 
cut on the negative real axis i.e. we integrate along s = e ±t7T x with x real 
and positive. Then one finds 



_d£_ 

2vri 



e s 



7-1 



sin(7T7) 



7T 



r(7) 



:i38i 



and so as z 



-00, 



6T(7) sin(7T7) A 



71 \ZP 



zP 



(139) 



The tail z — > +00: In this case, the asymptotic behavior of the integral 
in Eq. (|77j) can be evaluated by the saddle point method. The action hi(s) = 
—zs + 6s 7_1 inside the exponential has a saddle point on the positive real 
axis at s* = [z/bfy — l)] 1 ^ 7-2 . Since the function is analytic in the plane 
between the imaginary axes passing through s = and s = s*, the contour 
in Eq. (jTTj) can be shifted to s = s* and then the leading contribution to 
this integral for large z will come from the region around the saddle point at 
s = s*. Expanding hi(s) = hi(s*) + (s — s*) 2 h'[(s*) /2 + .. around the saddle 
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and performing the resulting Gaussian integration, one gets the positive tail 
of the scaling function as in Eq. (JHUj) . 



The value at z = 0: We now evaluate V^,(0) and show that its value is 
given as in Eqs. (fTH|) and (JEH). We need to perform the integral in Eq. (|77jl 
along the imaginary axis passing through s = 0. Rescaling s — > s6 -1 ^ 7-1 ^ we 
get 

/ lo ° A Q 
— e s7_1 (140) 
-too 2m 

To evaluate this integral we bend the contour along the rays s = exp^ 71 "^ 7-1 ) y 
(where y is real and positive) so that the argument of the exponential in ()140|) 
becomes real and negative. Then one finds 

ds_^-i = sin(7T/(7- 1)) I"™ d -a^r-i 
2lli 7T ./n 

r7^—l (in; 



7-1 V7 

where we have used J °° dxe~ xC " = T(l/a)/a and the identity pH4jl . Finally 
we have ^ 

VM = 6V(7-i) (7 _i)r((7-2)/(7-i))- (142) 
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